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Barnett Relaxation in Thermally-Rotating Grains 

A. Lazarian^ and W. G. Roberge^ 

ABSTRACT 

On ! We present an exact formulation of the physics of Barnett relaxation. Our 

formulation is based on a realistic kinetic model of the relaxation mechanism 
r^ ', which includes the alignment of the grain angular momentum in body coordinates 

^j ', by Barnett dissipation, disalignment by thermal fluctuations, and coupling of 

the angular momentum to the gas via gas damping. We solve the Fokker-Planck 
equation for the measure of internal alignment using numerical integration of 
the equivalent Langevin equation for Brownian rotation. The accuracy of our 
^ ! results is calibrated by comparing our numerical solutions with exact analytic 

cn '. results obtained for special cases. We describe an analytic approximation for 

_^ [ the measure of alignment which fits our numerical results for cases of practical 

O ! interest. 

ON 

Qh| Subject headings: dust, extinction — ISM, clouds — ISM, polarization 

O 

.. ' 1. Introduction 

>'■ 

^ ' The polarization of starlight, first discovered by Hiltner (1949) and Hall (1949), im- 

c^ • plies that the axes of dust grains are partially aligned with the interstellar magnetic field, B. 

Models of the Davis-Greenstein (DG) effect and other grain alignment mechanisms generally 
predict the alignment of the grain angular momentum, J, with respect to B (for a recent re- 
view, see Roberge 1996). However, to relate the predictions of these models to observations, 
one also requires the "internal distribution" of J in grain body coordinates. The classical 
papers on DG alignment (Jones & Spitzer 1967; Purcell & Spitzer 1971) assumed that the 
internal distribution is determined by gas damping and other external processes.0 However, 
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^For example, Jones & Spitzer (1967) noted that the component of B parallel to the grain angular velocity 
can enhance the alignment of J in grain coordinates. 



-2- 



Purcell (1979, hereafter P79) realized that frictional processes internal to the grain will dis- 
sipate rotational energy into heat on a timescale that is generally much shorter than the 
timescales for external interactions. Since they conserve angular momentum, these internal 
dissipation processes tend to drive the grain toward the state of minimum energy consis- 
tent with a fixed value of J, rotation about the principal axis of largest rotational inertia 
(henceforth the "major axis of inertia"). 

Several internal dissipation processes were discussed by P79, who found that "Barnett 
dissipation" is the dominant effect for typical interstellar grains. The mechanism is very 
elegant: Dolginov & Mytrophanov (1976) pointed out that a rotating grain has a magneti- 
zation M = — xi7/7 due to the Barnett effect, where x is the magnetic susceptibility, 7 is 
the magnetogyric ratio of the paramagnetic spins, and ft is the grain angular velocity. It 
is well known from theoretical mechanics (e.g.. Landau & Lifshitz 1976) that ft precesses 
in body coordinates unless the rotation axis coincides with a principal axis of the inertia 
tensor. Reasoning by analogy with ordinary paramagnetic dissipation in a time-dependent 
external magnetic field, P79 concluded that internal dissipation must also be caused by the 
time-dependent "Barnett equivalent magnetic field," -ffse = —^/l- 

Virtually all models of interstellar polarization developed since 1979 have assumed that 
Barnett dissipation aligns J perfectly with the major axis of inertia. However, Lazarian 
(1994, henceforth L94) has pointed out that this is a poor approximation if the grains rotate 
with thermal kinetic energies: thermal fiuctuations in the Barnett magnetization will excite 
rotation about all 3 of the body axes, preventing perfect alignment unless the rotation is 
suprathermal or the grain solid temperature is zero. 

The first model of internal relaxation to include thermal fiuctuations was presented by 
L94, who modeled the kinetics of Barnett relaxation by analogy with the DG effect. However, 
this analogy is not exact because, unlike Barnett relaxation, the DG effect involves external, 
nonconservative torques. This difficulty was ignored for the sake of simplicity by L94 and 
we now revisit the problem in order to carry out an exact analysis. In §2 we compare the 
efficiency of various dissipation mechanisms in order to delineate the circumstances under 
which the Barnett effect is the dominant process. In §3 we discuss Barnett alignment in an 
isolated grain whose angular momentum has a fixed but arbitrary magnitude, J. The value 
of J is determined by interactions of the grain with its environment and in §4 we formulate a 
kinetic model wherein the external interactions are provided by gas-grain collisions and the 
evaporation of molecules from the grain surface. We calculate the appropriate measure of 
internal alignment numerically in §5 and calibrate the accuracy of our results by considering 
special cases where exact solutions can be obtained analytically. Our results are discussed 
in §6 and summarized in §7. 
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2. Internal Dissipation 

Inelastic dissipation and Barnett dissipation were identified in P79 as the two major 
causes of internal dissipation. Inelastic dissipation is caused by the oscillating stresses asso- 
ciated with the motion of the grain angular velocity, CI, in body coordinates. It is possible to 
show that the timescale for inelastic dissipation, tai, scales as Q~^ compared to Q^"^ for the 
Barnett dissipation time, tear (see Eq. (|2-4|)). As a result, inelastic dissipation is the dom- 
inant process for sufficiently large angular velocities. However, for sufficiently rigid grains, 
these velocities lie in the range of suprathermal rotation; therefore, we consider only Barnett 
dissipation in the present paper, which deals with thermally rotating grains. The interesting 
problem of inelastic dissipation in dust grains will be discussed elsewhere. 

As we pointed out above, Barnett dissipation tends to align J with the major axis of 
inertia (see §1). Henceforth we will model the grains as oblate spheroids with semiaxes a 
parallel to z and b (with b > a) perpendicular to z, where z is the symmetry axis (=the 
major axis of inertia). Barnett dissipation decreases the angle 6 between J and z at a rate 

— = -G'Bar sine cos ^ (2-1) 

(P79; Roberge, DeGraff, & Flaherty 1993, henceforth RDGF93), where 

_ VKh\h-l)J^ 

^Bar - -^p , (2-2 j 

I z 

V is the grain volume, Ix is the inertia for rotation about an axis perpendicular to z, 
h = Iz/ 1x1 K is related to the imaginary part of the magnetic susceptibility at frequency Vt 
by 

.-.H^.Ox.0-.^(^). (.3) 

and Ts is the grain solid temperature. The scaling of K in equation (|2-3|) is appropriate for 
ordinary paramagnetic materials (see Draine 1996); for superparamagnetic grains, K should 
be increased by an uncertain factor ~ 10^. We define the Barnett time, tear, to be the value 
of Giar; SO that 

^""-^-^^ = VKhHh^W ^^'^^ 

The numerical value of tear depends on J; a typical angular momentum for thermal rotation 
at the gas temperature, Tg, is 



and the corresponding numerical value of tear is 

tBa.( Jtw.) = 4 X 10^3^0 T,: KZI, b\ ^"^(|^_l;^,J^f^^'^' s, (2-6) 
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for a grain with radius 10^^6_5 cm and density pso g cm"^. We may say that eq. 
describes Barnett dissipation for an individual grain with fixed angular momentum J, while 
eq. ( ^-(j| ) decsribes Barnett dissipation for an ensemble of grains with rotational temperature 



3. Barnett Dissipation in an Isolated Grain 

Barnett dissipation is able to transform rotational energy into heat because coupling ex- 
ists between the grain rotational and vibrational degrees of freedom. However, the existence 
of this coupling implies that energy can also be transferred in the opposite direction, i.e., 
from vibration to rotation (see Jones & Spitzer 1967). As a result, fiuctuations in the thermal 
energy will randomly excite rotation about all three body axes with energies of order fcTg. If 
the rotational energy of the grain is ^ kTg, then these thermal fiuctuations are negligible^ 
and Barnett dissipation aligns J perfectly with z. This is the assumption adopted by P79, 
who was concerned only with rotation at highly superthermal energies. However, the grains 
in molecular clouds are likely to have rotational energies ~ kTg. Consequently, it is of great 
practical interest to determine the distribution of orientations established by the combined 
effects of Barnett dissipation and thermal fluctuations. 

This distribution can be calculated analytically for an isolated grain. Let /te (^ | J) be 
the distribution function for 6, deflned so that /te sinO d6 is the probability of flnding J 
oriented between 6 and 6 -\- dO. Here we regard /te to be a function of 6 with J a flxed but 
arbitrary parameter; the determination of J by external interactions is discussed in §4. Since 
the grain is assumed to be an isolated system in thermodynamic equilubrium, /te can be 
found, in principle, by maximizing the entropy of the combined rotation/vibration system. 
However, the calculation is drastically simplifled by noting that the number of vibrational 
degrees of freedom (~A^, where A^ is the number of atoms in the grain) is vastly greater than 
the number of rotational degrees of freedom (=3). Thus, the rotational degrees of freedom 
may be regarded as a small system immersed in a heat bath at constant temperature, Tg.^ 



^However, see Lazarian &; Draine 1997, where it is shown that the effects of incomplete Barnett ahgnment 
are important during "crossovers" in suprathcrmally rotating grains. 

•^Here we assume for simpUcity that Tg is maintained at a constant value by the various processes (inter- 
actions with the ambient radiation field, gas damping, etc) by which energy is exchanged between the grain 
and its environment. That is, we assume that the radiation field, gas temperature, etc change on a timescale 
that is much longer than tear- 
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It follows that /te is just the Boltzmann distribution, 

/te (^ I J) = C exp[-^,ot W ATs] 
where C is a normalization constant and 



ErotiO) 



2L, 



1 + (/i - 1) sin^ e 



(3-1) 



(3-2) 



is the rotational energy. 



In a steady state, the linear and circular polarization produced by an ensemble of 
spheroidal grains depend only on the measure of internal alignment. 
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(3-3) 



where angle brackets denote the distribution average. Evaluating the average for /te, we 
find that the measure of internal alignment for an isolated grain is 

/o« t' exp (t2 - e) dt 



Qx,TE{a/b,Ts,J) 



2e 



J^ expit^-e)dt 



(3-4) 



where we have introduced the dimensionless "suprathermality parameter," 

The dependence of Qx,te on ^ (Fig. 1) has a straightforward physical interpretation. For 
a fixed value of J, the difference in rotational energy between the minimum {6 = 0) and 
maximum {6 = tt/2) rotational energy "states" of a grain is AE = {h — 1) 7^/2/^ [see eq. 
(P^)]. Thus, ^^ is just the ratio of the energy required to disorient J to the energy added 
to rotation by a typical thermal fluctuation. Weak alignment (^ <^ 1) occurs when a typical 
fluctuation is large enough to disorient J [kTs ^ AE), corresponding to nearly-spherical 
grains {h ^ 1) or thermal rotation energies (J^/2/2 ~ kTg). Strong alignment {^ ^ 1) 
occurs only for nonspherical grains and suprathermal rotation. 



4. Barnett Dissipation with External Interactions 
4.1. Mathematical Description 

The effects of external interactions can be included by solving the Fokker-Planck equa- 



tion. 



dt 



-V-5, 



(4-1) 
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where t is time, / (J, 9) is the joint distribution for J and 9, and S is the probabihty current. 
The latter is given in general by 

S = Af - l-V-{Bf), (4-2) 



where 



is the mean torque and 



2 

is the diffusion tensor. In principle, A and B include the cumulative effects of every process 
that changes J in the body frame. In this paper, we consider the effects of Barnett dissipation 
and gas damping but neglect the effects of other external interactions such as the Davis- 
Greenstein mechanism. As a result, our calculations accurately describe the internal angular 
momentum distribution only under conditions where the gas damping time is much smaller 
than the timescale for other external processes.^ We assume that the grain center of mass is 
at rest with respect to the gas; the generalization of the results presented here to Gold-type 
alignment will be considered in subsequent papers (Lazarian 1997b; Roberge & Lazarian 
1997b). We also assume that the properties of the grains and their environment are steady 
on timescales much longer than the gas damping and Barnett times, so that only the steady 
solution of the Fokker-Planck equation is of interest. 



4.2. Diffusion CoefRcients for Barnett Relaxation 

It is simplest to calculate A^^^"^' and B^^^^\ the mean torque and diffusion tensor asso- 
ciated with Barnett relaxation, in a spherical polar coordinate system with basis (j, 0, 0j, 
where the polar axis is parallel to the grain symmetry axis. The conservation of J by the 
Barnett effect and the azimuthal symmetry of an oblate spheroid together imply that the 
only nonvanishing components are Ag '^'^ and Bqq^"^ . The former follows immediately from 
equation ( |2-1|) and is given by 

^(Bar) _ ^ sin e COS e (4-5) 

^Bar 



''in particular, we do not attempt to describe magnetic dissipation by ordinary or superparamagnetic 
grains. These effects will be considered in subsequent papers on Davis-Greenstein alignment (Lazarian 
1997a; Roberge & Lazarian 1997a), which generalize the treatments of DG alignment given in Roberge et 
al. (1993) and Lazarian (1995a). 
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where tear is defined by eq. ( P-4|) . The analogous expression for Bqq^^ can be determined 
by noting tliat, in tfie absence of otfier processes, Barnett relaxation must establish the 
thermodynamic equilibrium distribution [eq. ( |3-1| )]. Further, the rate of every microscopic 
process equals the rate of the time-reversed process in thermodynamic equilibrium (the 
"principle of detailed balance"), so the probability current must vanish at every point in 
phase space. Writing out the components of S in polar coordinates and setting Sq = 0, we 
find that detailed balance requires 



sin e /te ^r'''^ ~ 2Jde [^^'^ ^ •^'^^ ^' 



.(Bar) 



(Bar) 



0. 



Expression ( [4-6| ) is a first-order differential equation for Bqq^^ . The solution is 



B, 



(Bar) 



J2 



exp(^^sin^^) / v^exp(-^\)d|/ + exp(-^ 



tear sin e 
In solving equation (|4-6|) , we have imposed the condition^ 

p(Bar),7r _ J^ 

which insures that the Bg^^'^ is smooth at '7r/2: 



2 2 

COS I 



(4-6) 



(4-7) 



(4-8) 



(4-9) 



When J^ <^ IzkTg thermal fluctuations easily change its direction and this causes Bqq^^' -^ cx) 
in agreement with our results. 



4.3. Diffusion CoefRcients for Gas Damping 

The diffusion coefficients for gas damping were calculated by RDGF93 for oblate spheroids 
on the assumption that f2 is aligned perfectly with z (i.e., perfect Barnett alignment). How- 
ever, the diffusion tensor is independent of ft to first order in Qb/vth_, where Vtt ■ 



'2kTg/m 

is the gas thermal speed. We may therefore adopt the results of RDGF93 for the diffusion 
tensor with high accuracy. The nonzero Cartesian components are 



BiT^ = ^-^nmb^lT^ie) 



Z 

n 



(4-10) 



^It worth noting that condition (4-8) coincides with the result for the diffusion coefficient at a critical 
point obtained in Zeldovich (1942) for the problem of phase transitions of the first order. 



and 



D(gas) _ o(gas) 

yy XX ' 



5lr^ = ^^^&\\r||(e)[i + ^' 



(4-11) 



(4-12) 



where n is the gas number density and the geometrical factors rn and F^ are dimensionless 



(4-13) 



functions of a/h. They are given in terms of the eccentricity, e = Jl — a^/fe^^ by 

3 



rii(e) = tt: { 3 + 4(1 - ^)g{e) - e"^ [l - (1 - ^fg{e)\ } 



and 



r^(e) = A{7-e^ + (l-e2)^^(e) + (l-2e2 



32 
respectively, where 



1 + e 



-2 



qie) = — In ( 

^^ ^ 2e VI -e 



l-(l-e2)Ve)]]}, (4-14) 



(4-15) 



The mean torque due to gas damping vanishes to zeroth order in 17, so we cannot adopt 
the expressions given by RDGF93. However, one can relate the mean torque to the diffusion 
tensor via the principle of detailed balance^ with the result 



and 






^(gas) ^ 4v^^™^4 



3L 



nmh VtiXwJz- 



(4-16) 
(4-17) 

(4-18) 



If the Barnett alignment is perfect, then J^ = Jy = and the preceding expressions reduce 
to the components obtained in RDGF93. The characteristic timescale for gas damping is 



t 



SL 



gas 



A^/^Tn'mb^VtlJP\\ 



where 

for a gas of pure H2. 



igas = 7x10^ pso ^4 ^ T 1^ a-5 Ty s 



(4-19) 
(4-20) 



^Strictly speaking, this principle applies only in thermodynamic equilibrium with Tg = Tg. However, it 
can be shown that the mean torque is independent of Ts/Tg to first order in flb/vth so that, apart from a 
small error, the thermodynamic equilibrium result equals the result for Ts ^ Tg. 
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4.4. Numerical Methods 

We solve the Fokker-Planck equation indirectly by integrating the equivalent "Langevin 
equation for Brownian rotation," 

dJ = Adt + y/BdW (4-21) 

(RDGF93), where ^/B denotes the matrix square root of B and dW is a vector of statistically 
independent, Gaussian random variables with variance dt. To calculate Qx, we integrate 
equation ( |4-21D numerically to generate a simulation of J{t) and equate the ensemble average 
of cos^ 6 in eq. (|3-3|) to the time average over the history of the simulation. The numerical 
integration was performed using the Euler algorithm; this scheme approximates the change 
in Ji during each time step. At, by 

AJi = Ai At + (/b) . . AWj, (4-22) 

where the AWj are independent samples from a Gaussian distribution with variance (At) . 
The numerical integration was performed using a "mixed" scheme, where the changes in 
Jx, Jy, and Jz due to gas damping were computed from the Euler algorithm in Cartesian 
components and the change in 6 due to Barnett relaxation was computed from the Euler 
algorithm in polar coordinates. 



5. Results 

5.1. Dimensionless Parameters 

The number of independent parameters in the problem is greatly reduced by rewriting 
the Fokker-Planck equation in terms of dimensionless time and angular momentum variables. 
With time measured in units of tear and angular momentum measured in units of Jtherm, 
only 3 dimensionless parameters appear: 

a/b = grain axis ratio, (5-1) 

^ ,^ dust solid temperature 

Ts/T, = , (5-2) 

gas kmetic temperature 

and 

Sb - ^ 'T y (5-3) 

''Barl^'-'thermJ 

Henceforth we will refer to 6b as the "Barnett damping parameter." In the model adopted 
here, where gas damping is the only external interaction, Qx is completely determined by 
specifying the values of a/b, Tg/Tg, and 6b- 
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5.2. Benchmark Tests 

The accuracy of our numerical method can be calibrated by considering two special 
cases where exact analytic solutions for Qx exist. Consider first the case of an isolated 
grain, for which Qx = Qx.te- To simulate an isolated grain, we assigned J/Jtherm as an 
initial condition, turned off gas damping in our numerical code, and integrated the Langevin 
equation with the timestep and averaging time set to At = 1 x 10^^ and T = 2 x 10^ 
dimensionless units, respectively. The measures of alignment computed for several cases 
are plotted as solid circles in Figure 1. The rms error in Qx is 4 x 10~^ for the numerical 
solutions in Figure 1. 

Exact solutions can also be obtained in the "Maxwellian limit," 5-q = 0, where the 
timescale for Barnett relaxation is infinitely long compared to the timescale for gas damping. 
The Maxwellian distribution for 6 is 

fuUO) = ^ (cos^ 9 + h sin^ ey'^' (5-4) 

(Jones & Spitzer 1967, but note that we have corrected a minor typographical error) and 
the corresponding measure of alignment is 



Q 



X,Max 



2il-h-^] 



1 ; sin"^ (l - h- 



v^ 



\ (5-5) 



In the Maxwellian regime, the partial alignment of J with z is due entirely to the "inertial 
asphericity" of the grains and, appropriately, Qx depends only on a/6.[] In Figure 2, we 
compare the exact solution for Qx.Max (solid curve) with numerical solutions calculated for 
a few values of ajh (filled circles). The numerical solutions were obtained by integrating the 
Langevin equation as a function of time in units of tgas (since tear is irrelevant for this case), 
with At = 1 X 10^^ and T = 2 x 10^ dimensionless units, respectively. The rms error in Qx 
for the filled circles in Figure 2 is 2 x 10^'^. 



5.3. Barnett Relaxation with Gas Damping 

Typical results for realistic calculations including gas damping and Barnett relaxation 
are shown in Figures 3 and 4. As noted above, the case 5^ = corresponds to the Maxwellian 
regime and, appropriately, all of the solutions in Figures 3 and 4 converge to the Maxwellian 
case as bs —>■ 0. As 6b increases, Barnett relaxation becomes increasingly important and 



^Recall that h — 2/ (l + a^ /b^) for a homogeneous, oblate spheroid. 
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Qx diverges from Qx.Max- Notice that Barnett relaxation always enhances the alignment, 
in the sense that Qx > Qx,Max for 6b > 0. For large Sb, the numerical solution for each 
combination of a/b and Tg/Tg approaches a limiting value 

QxM l*^/^' TjTg) = lim Qx (a/b, TjTg, 5b) , (5-6) 

which depends only on the grain shape and dust-to-gas temperature ratio. This limit rep- 
resents the maximum possible enhancement of the alignment (relative to Qx,Max) which can 
be produced by Barnett relaxation in a thermally-rotating grain. The qualitative behavior 
of the results presented in Figures 3 and 4 is consistent with the predictions of L94. For 
example, the maximum efficiency of Barnett alignment decreases as Tg/Tg increases and as 
the grains become more spherical. 

Note that the oscillations which occur for large 6b in Figs. 3 and 4 are spurious, that 
is, they are too large to be statistical fluctuations. The oscillations are probably an artifact 
of the ad hoc prescription we use to treat the singularities in Bgg^^ at ^ = and 9 = n [see 
eq. (p^)]: in computing Bg^^^' , we set 6 = 6'min whenever the trajectory "visits" the interval 
(0, ^min) and similarly for points in the interval (tt — ^min, tt). This procedure is motivated by 



the fact that the fluctuations in 9 are proportional to y Bqq^^ [see eq. ( [4-2 1|) ]. For some value 
of ^min, the fluctuations for < O^^^^ become so large that the grain is completely disoriented 
and, roughly speaking, the precise values of the fluctuations ought to be irrelevant. After 
experimenting with different choices, we set ^min = 0.1 rad for the calculations presented in 
Figs. 3 and 4. 

The fact that the solutions in Figures 3 and 4 approach limiting values for large 5b is 
a consequence of elementary thermodynamic considerations. The case 5^ 3> 1 corresponds 
to conditions where 6 evolves rapidly (on a timescale ~ tear) due to the effects of Barnett 
relaxation while J evolves slowly (on a timescale tgas 3> tear) due to gas damping. The 
disparity between tear and tgas means that, insofar as the motion of is concerned, we may 
regard J as a slowly-varying parameter. Thus, if we compute the measure of alignment 
by averaging the motion of 9 over times much longer than tear but much shorter than 
tgas, the result will be Qx ~ Qx,te, the approximation becoming exact in the limit 5b -^ 
cxD. We conclude that the asymptotic behavior of the solutions in Figures 3 and 4 occurs 
because, when 5b 3> 1, the distribution of 6' is always close to the thermodynamic equilibrium 
distribution and the latter is independent of 5^. Accordingly, we will refer to the case 5b ^ ^ 
as the "thermodynamic regime." Of course, Qx,te depends on the value of J, so the value of 
Qx computed by the "short-time averaging" procedure described above will change quasi- 
statically as gas damping slowly changes the magnitude of the angular momentum. Thus, it 
is necessary to compute Qx by averaging the motion over many gas damping times, as we 
have done for the calculations shown in Figures 3 and 4. 
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For real interstellar grains, the Barnett time is typically a few orders of magnitude 
shorter than the gas damping time [compare eqs. ( p^ ) and ( [4-2(]| )] so the solutions of practical 
interest in future applications correspond to 6b ~ 10^-10'^. Although we are unable to 
compute solutions for 5b values this large (because the Langevin equation becomes too 
stiff), Figures 3 and 4 suggest that the solutions are already close to the limit (5b ^ cxd when 
5b = 25. Accordingly, values of Qx computed for 6b = 25 should be good approximations 
to the results of practical interest. We have calculated the measure of internal alignment 
as a function of a/b and Ts/Tg for the case 6b = 25; the resulting values of Qx appear as 
the upper entries in Table 1. The entries for Ts/Tg = 1 provide an additional check on our 
calculations because the angular momentum distribution must be Maxwellian when all of 
the temperatures in the system are equal. We have plotted these entries as open circles in 
Figure 2; the rms error in Qx is 6 x 10"'^ for these points. 

Because the distribution of 6 is always close to the thermodynamic equilibrium distri- 
bution when 6b is large, one should be able to approximate the values of Qx in Table 1 by 
an expression of the form 

Qx~Qx,TE (a/6, T„ Jeff), (5-7) 

where Jes is some "effective J-value." We have computed the effective J value for each entry 
in Table 1 by interpreting expression (|5-7| ) to be the definition of Jeff- In particular, each 
Qx value in the table is given exactly by 

Qx = Qx,TE {a/b, Ts, a Jiviax) , (5-8) 

where 



J 



Max 



^ , ja/by 



Ts 

1 + 777 



1/2 

'-'therm v"-*"^/ 



is the rms value of J for a Maxwellian angular momentum distribution and a is a dimen- 
sionless number. For each parameter combination, the lower entry in Table 1 gives the 
corresponding value of a. Figure 5 gives a graphical comparison between selected results 
from the table and the "a = 1" approximation The agreement is generally good, although 
the a = 1 approximation systematically underestimates Qx for points with Qx ^ 0.4. This 
systematic discrepancy is likely due to errors in the numerical solutions associated with the 
singularity in Bgg^'^ at 9 = rather than the quality of the approximation. In view of the 
various uncertainties in grain properties, e.g., in the rotational temperature and rotational 
inertias, we regard the a = 1 approximation as a reasonable method for estimating the degree 
of Barnett alignment in real grains when gas damping is the dominant external interaction. 
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6. Discussion 

The study of internal relaxation presented in this paper includes the physics of Barnett 
relaxation and gas damping but does not include the Davis- Greenstein effect or Gold's mech- 
anism. A quantitative study including the omitted processes is a challenging problem, which 
we address in subsequent papers in this series (Lazarian 1997a and Roberge & Lazarian 1997a 
for the Davis-Greenstein mechanism and Lazarian 1997b and Roberge & Lazarian 1997b for 
Gold's mechanism). Without going into details, we may summarize the qualitative effects 
of incomplete Barnett alignment on the DG and Gold mechanisms as follows. Because the 
DG and Gold mechanisms are both concerned with the alignment of J with respect to B, 
the wobbling of a thermally-rotating grain about its angular momentum (which occurs when 
Barnett alignment is imperfect) will result in a decrease in its polarizing efficiency relative 
to the case for perfect Barnett alignment. According to calculations based on the assump- 
tion of perfect Barnett alignment (Lazarian 1994; Roberge, Hanany, & Messinger 1995), the 
efficiency of Gold alignment depends on the orientation dependence of a grain's cross section 
for collisions with streaming gas particles, in the sense that the alignment becomes more 
efficient as this difference increases. The precession of the grain axes about J which occurs 
when Barnett alignment is imperfect will generally reduce this difference. Consequently, 
we anticipate that imperfect Barnett alignment will reduce the efficiency of DG and Gold 
alignment for thermally-rotating grains. Conversely, the DG and Gold mechanisms affect 
the efficiency of Barnett alignment. Paramagnetic dissipation in the interstellar magnetic 
field will decrease the Barnett alignment by reducing the value of J . On the other hand. 
Gold's mechanism increases J and therefore increases the efficiency of Barnett alignment. 
Strangely, incomplete Barnett relaxation increases the efficiency of paramagnetic alignment 
if the grains rotate suprathermally by increasing the correlation in the directions of J before 
and after crossover events (Lazarian & Draine 1997). 



7. Summary 

The principal results of this paper are as follows: 

1. We have presented an exact formulation of the physics of Barnett relaxation which 
includes the effects of thermal fluctuations in the grain material. 

2. We have shown that the measure of internal alignment, Qx-, depends only on 3 di- 
mensionless parameters in the approximation where the external torques on the grain 
are dominated by gas damping. These parameters are the grain axis ratio, a/6, the 
dust-to-gas temperature ratio, Tg/Tg, and a "Barnett damping parameter," 5^. 
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3. We have calculated Qx for a large number of parameter combinations. 

4. We have obtained an analytic approximation that accurately reproduces our numerical 
solutions in the portion of parameter space relevant to studies of interstellar polariza- 
tion. 

5. The results of this paper confirm the conclusion of L94, that the alignment of J with 
respect to the axis of major inertia is incomplete for thermally-rotating grains. The 
incomplete alignment stems from thermal fluctuations within the grain material and 
becomes increasingly significant as the grain temperature rises. The rotation of a grain 
with J parallel to its major axis of inertia corresponds to the minimum kinetic energy 
for a given J value. Thermal fluctuations disorient J because they transfer energy 
from the grain solid material to the rotational degrees of freedom. In a steady state, 
energy is transferred back to the grain material at an equal rate by Barnett dissipation. 

6. It is important to include incomplete Barnett alignment in studies which relate grain 
alignment theory to polarimetric data. The polarization of electromagnetic radiation 
is caused by the preferential alignment of the grain axes. At the same time, theoretical 
studies of alignment deal with anisotropics in the angular momentum distribution.^ 
The distribution of J in grain coordinates is essential for relating observations to 
theory. 
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FIGURE CAPTIONS 

Fig. 1 — Measure of alignment for an isolated grain with constant angular momentum, 
plotted versus the square root of the suprathermality parameter. Solid curve: exact solu- 



tion from equation (|3-4|) . Filled circles: approximate solutions obtained by integrating the 
Langevin equation. The rms error in the numerical solutions is 4 x 10~^. 

Fig. 2 — Measure of alignment in the Maxwellian regime plotted versus the grain axis ratio. 
Solid curve: exact solution from equation (p-5|). Filled circles: numerical results for 6b = 
and arbitrary Tg/Tg. Open circles: numerical results for 6b = 25 and Tg/Tg = 1. The rms 
error in the numerical solutions is 2 x 10~^ for the filled circles and 6 x 10~^ for the open 
circles. 

Fig. 3 — Numerical calculations of the measure of alignment for a highly-flattened grain 
with a/b = 0.1 including the effects of gas damping and Barnett relaxation. Abcissa: Bar- 
nett damping parameter, 6b [see eq. (^^)]. Ordinate: Measure of internal alignment, Qx- 
Solutions are presented for Tg/Tg = 0.1 (filled circles) and Tg/Tg = 0.5 (open circles). The 
dashed line is the exact solution in the Maxwellian regime and corresponds to Tg/Tg = 1. 

Fig. 4 — Similar to Figure 3 but for a modestly flattened grain with a/b = 0.5. 

Fig. 5 — Graphical comparison between the analytic approximation defined by expression 
( p-8| ) with a = 1 (solid curves) and selected numerical results from Table 1 (symbols). 
Results are shown for a/b = 0.1 (top curve, filled circles), a/b = 0.5 (middle curve, open 
circles), and a/b = 0.9 (bottom curve, filled squares). 
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Table 1 
Measure of Internal Alignment for Sb = 25 













Ts/T, 










a/b 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


0.1 


0.54 


0.39 


0.29 


0.24 


0.22 


0.18 


0.18 


0.16 


0.14 


0.13 




0.85 


0.95 


0.96 


0.98 


1.01 


0.98 


1.03 


1.00 


0.98 


0.97 


0.3 


0.49 


0.35 


0.25 


0.22 


0.18 


0.17 


0.15 


0.14 


0.14 


0.13 




0.85 


0.95 


0.95 


0.98 


0.98 


0.99 


0.99 


1.00 


1.02 


1.01 


0.5 


0.44 


0.28 


0.21 


0.18 


0.15 


0.13 


0.12 


0.11 


0.10 


0.09 




0.90 


0.97 


0.98 


1.04 


1.02 


1.00 


1.00 


1.01 


1.00 


0.96 


0.7 


0.30 


0.17 


0.13 


0.11 


0.09 


0.08 


0.07 


0.07 


0.06 


0.06 




0.93 


0.96 


1.01 


1.02 


1.00 


0.98 


0.97 


1.01 


0.99 


0.98 


0.9 


0.11 


0.06 


0.04 


0.03 


0.03 


0.02 


0.02 


0.02 


0.01 


0.01 




0.98 


0.97 


0.97 


0.97 


0.94 


0.93 


0.88 


0.84 


0.82 
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